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ABSTRACT 

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 
GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. 
From a piecewise Gaussian ht to each flare, time lags At between the observation frequencies 
1 / and spectral indices a based on peak amplitudes A are determined. From the fit At cx 
= 1.10 ± 0.18 indicating equipartition between the magnetic field energy density 
and the particle energy density. From the fit A c>c a is in the range —0.24 to 1.52. A 
mean magnetic field strength at 1 pc, Bi = 0.5 ± 0.2 G, and at the core, Score = 46 ± 16 
mG, are inferred, consistent with previous estimates. The measure of core position offset is 
flri, = 6.4 ± 2.8 pc GHz^/^’’ when averaged over all frequency pairs. Based on the statistical 
trend shown by the measured core radius Vcore as a function of z/, we infer that the synchrotron 
opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law 
slopes in the range —1.6 to —3.5 describing the power spectral density shape and gives bend 
timescales in the range 0.52 — 0.66 yr. This result, and both positive and negative a, indicate 
that the flares originate from multiple shocks in a small region. Important objectives met in 
our study include: the demonstration of the computational efficiency and statistical basis of 
the piecewise Gaussian fit; consistency with previously reported results; evidence for the core 
shift dependence on observation frequency and its utility in jet diagnostics in the region close 
to the resolving limit of very long baseline interferometry observations. 


Key words: galaxies: active - (galaxies:) quasars: general - (galaxies:) quasars: individual: 
3C 454.3 - methods: data analysis 


1 INTRODUCTION 

Two groups of objects, BL Lacertae (BL Lac; with featureless op¬ 
tical spectra) and flat spectrum radio quasar (FSRQ; having promi¬ 
nent emission lines) constitute a subset of radio-loud Active Galac¬ 
tic Nuclei (RL-AGNs) called blazars, which according to the ori¬ 
entation based unified schemes of RL AGNs have jets making an 
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angle of ^ 10° or so to the line of sight (e.g. |Urry & Padovani] 
|1995| >. AGNs are characterized by strong flux variability at all 
wavelengths, strong and variable polarization in the radio to optical 
wavelengths, and a core dominated radio spechum. Blazar emis¬ 
sion from the relativistic jet is dominated by non thermal emis¬ 
sion (synchrotron and inverse Compton) which is Doppler boosted 
making a cone angle of ~ l/F where F is the bulk Lorentz fac¬ 
tor. Flicker in radio flux from these AGN can be caused by inter¬ 
stellar scintillation. The apparent extremely high brightness tem¬ 
peratures of > lO^^K inferred from radio observations can be 
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attributed to relativistic beaming, coherent radiation mechanisms 
and geometric effects jWagner & Witzel|1995^ . Periodic variabil¬ 
ity in the blazar might he caused hy the helical movement of the 
jet components, shocks, jet precession effects, accretion disc in- 
stahilities jHonma et al.|19 91[|Camen zind & Krockenberger 1992 
Gopal-Krishna & Wiita|1992[|G6mez et al.|1997[|Bach et al. 2006 


Bri tzen et al.|2 010 1 or due to the rotation of the secondary hlack 

hole around the primary supermassive black hole (e.g. |Sillanpaa| 
|et al.|p^ 88 ] [Valtonen et al.|[2008] [Rieger & Mannheim||2000| ^ 
Doppler boosting shortens the variability timescales by a factor of 
5“^ where this Doppler factor depends on the Lorentz factor as 
5 = l/r(l — PcosO) where 0 is the angle between the jet direc¬ 
tion and the observer line of sight and /? = (! — is the 

speed of the jet component in the source frame. 

3C 454.3 (PKS 2251-1-158) is a highly variable blazar, classi- 
fied as a FSRQ with red shift z = 0.859 (e.g. [Hewitt & Burbidgej 
with a black hole mass of 4.4 x 10^ Mq 


1989 


Jorstad et al. 2001 


I Gu, Gao & Jiang|2001 1 . An ultra-violet excess due to the domi¬ 


nance of thermal emission from the accretion disc over synchrotron 
radiation is inferred in jRaiteri et al.|j2b07||2008) . Since 2001, the 
FSRQ has displayed strong variability over the entire electromag¬ 
netic spectrum. It underwent a bright outburst almost simultane¬ 
ously from radio to X-ray bands in May 2005 when it reached its 


2006 

Vercel- 

2008] 

Jorstad 


jet al.| 2010 ^ . Flux variations at 230 GHz were detected around two 
months after the optical variations. This outburst was inferred to 
be caused by shocks propagating along the relativistic jet when 
the emitting region at optical wavelength passes closest to the ob¬ 
server’s line of sight first and does so later on for other wavelengths 
( jVillata et al.|2007^ . The 2005 outburst was followed by a compar¬ 
atively quiescent period till spring 2007 at all frequencies, when 
it displayed a “big blue bump”, characteristic of radio-quiet AGNs 
and indicating thermal emission from the accretion disc ( jRaiterij 
jet al.|2007l l. |Katarzyhski & Ghisellinl| ( |2007^ proposed a model to 
explain the activities of the FSRQ during 2005-2007 time period, 
according to which the location from which bulk kinetic energy 
is dissipated along the jet is dependent on the bulk Lorentz factor 
and compactness of the perturbations which arise. Correlations be¬ 
tween light curves at different frequencies during the 2005-2006 
flare established that a single radiation mechanism is responsible 
for both radio and optical emission. Observations were carried out 
on all possible timescales, i.e., from days to years, in both optical 
and radio from which it was found that the duration of the flare and 
the individual features of the flare were same in both, thus showing 
that both were produced by the Doppler boosted relativistic jet. 

A multiwavelength monitoring campaign during a high flux 
state in 2008 jBonning et al.|[2009^ indicates a strong correlation 
between infrared, optical, ultraviolet and 7 -ray wavelengths. This 
is attributed to the external Compton scattering where accretion 
disk and emission line photons are scattered by strongly relativis¬ 
tic electrons (with Lorentz factors 10^ — 10^) which are also ra¬ 
diating synchrotron based infrared and optical emission. A study 
of intra-day variability in the source during 2009 - 2010 jGaur,| 
[Gupta & Wiita|[2012| ) during a high flux state indicates a strong 
correlation between the 7 -ray and optical wavelengths and that the 
former leads the latter by 4.5 ±1.0 days and is attributed to exter¬ 
nal Compton scattering. A multi frequency study from millimeter 
wavelengths to 7 -rays during three outburst events in 2009 - 2010 
( [Jorstad et al.|201^ reveals a similar radiation mechanism and lo¬ 
cation of the flares with small to no time lag being found between 
optical and 7 -ray light curves. Typical variability timescales are in¬ 


ferred to be ~ 3 hours in the parsec scale jet adding to the conclu¬ 
sion that the emission is highly localized. The measured magnetic 
field strength is in the range 0.1 — 1 G. A study to derive the bulk 
kinetic power using jet parameters derived in the framework of the 
[Konigl[ jl981| l inhomogeneous jet model at the 43 GHz observa¬ 
tional frequency indicates a jet viewing angle of 3.8°, bulk Lorentz 
factor of 32.4 and Doppler factor of 20.3 jGu, Cao & Jiang|2009| (. 
From the derived bulk kinetic luminosity and the luminosity from 
the broad line emission (due to disk based sources), the study in¬ 
fers a strong correlation between these two indicating that the jet 
power is intrinsically linked to the accretion process. Owing to its 
high flux state and broadband spectral energy distribution (SED), it 
offers opportunities for intra-day variability (IDV) studies ( |Tavec-[ 
[chip et al.|2010| ) as well as spectral variability studies. A study of 
the relationship between the bulk Lorentz factor and the AGN prop¬ 
erties ( |Chai, Cao & Gu|2012^ estimates jet parameters derived in 
the framework of the inhomogeneous jet model ( |Konigl|1981[ l. A 
magnetic field strength at 1 pc of 0.493 G was obtained, and from 
the large sample of radio-loud AGN, a strong correlation is inferred 
between the bulk Lorentz factor and the black hole mass, indicating 
that accretion and the resulting black hole spin up could play a role 
in jet power. 

In section|^ we describe parameters that can be derived from 
the analysis of the radio light curves, including the magnetic field 
strength and size of emitting core, based on physical models pro¬ 
posed in earlier literature. In section]^ we describe the observations 
made to obtain the final analyzable light curves and the analysis 
procedure which includes a piecewise Gaussian fit to each flaring 
portion in the light curve. In this section, we then describe the time 
series analysis of each light curve segment using the Fourier peri- 
odogram. In section|^ the results of the analysis and the estimated 
parameters including the time lags, spectral indices, magnetic field 
strength, size of emitting core and others are presented. This is fol¬ 
lowed by a discussion of the results from the time series analysis. 
We finally conclude in sectionj^with the comparison of our results 
with earlier work and its impact. 


2 FREQUENCY-DEPENDENT CORE SHIFTS 


Due to finite VLBI resolution, the true optically thick radio core 
(i.e. surface with optical depth r = 1 ) is not completely resolved. 
Since the t — 1 surface is frequency dependent ( |Konigl[ri981[ 
[O’Sullivan & Gabuzda|200^ , the radio cor e position dep ends on 
the observation frequency t/ as r oc iKonigl 1981 1 , with r 

being the distance from the central region and the index hr given 
by jLobanov|1998| l, 


kr = 


(3 — 2a)m ± 2n — 2 


( 1 ) 


5 - 2q; 

under the assumption that the ambient medium pressure on the 
jet is negligible ( |Konigl|1981[ |Lobanov[[ 1 9^ and external pres¬ 
sure is non-negligible with non-zero gradient along the jet. In the 
above expression for a is the spectral index (S oc iz“) while 
m and n describe the power law decay of the magnetic field (B — 
i3i(r/lpc)“™) and particle number density (N = A^i(r/lpc)“") 
respectively, with distance from the jet base (i.e. the jet apex). Here 
Ni and Bi refer to values of N and B at a distance of r = 1 
pc from the core in the observer’s frame. These power law expres¬ 
sions are valid till ~ 1 pc as the jet shape changes from conical 


to parabolic closer to the true jet origin (e.g., Krichbaum et al. 


[20061 [Asada & Nakamura||2012( [Nakamura & Asada||2013(^ 

[skin & Nokhrina|2006[ [McKinney|2006^. Assuming equipartition 
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between electron energy density and the magnetic energy density 
we obtain n = 2m ( |Lobanot^|1998) . Assuming the combination 
m = 1 and n = 2 which describes well the synchrotron emission 
from compact VLBI cores ( |Konigl| 1981 [ |Lobanov| 1998) > for con¬ 
stant jet flow speed along with constant opening angle, we obtain 
fcr = 1, independent of the spectral index a. 

In the shock-in-jet model scenario, shocks are produced as a 
result of change in physical quantities such as velocity, electron 
density, magnetic field through the jet or pressure at the base of 
the jet. These shocks propagate down the jet and cause core bright¬ 
ening, thus resulting in a flaring state of the source which is fol¬ 
lowed by the appearance of a jet component in VLBI images (e.g. 
[Marscher & Gear|I985[ [Gomez et al.|1997| l. Let Ro„ be the dis¬ 
tance of the onset of a shock which crosses the r = 1 surface at 
time Ti for frequency Vi thus giving a flare in the light curve. The 
core position shift with frequency gives time lags between differ¬ 
ent frequency light curves. For apparent superluminal motion (e.g. 
|Rees| 19671 [Tiirler, Courvoisier & Paltani|2000^ , the time t in the 
observer’s frame passed after the onset of the shock wave at the 
distance Ron is ( [Kudryavtseva et al.|2011| ) 

^ Ro^{u)) 

^appC 

where R is the distance along the jet axis, ftapp is the apparent 
superluminal velocity of the jet component and 6 is the jet viewing 
angle. 

The time lag between light curves at two frequencies is given 
by 


At = t{Va) — t{Vb) 

_ (l-|-z)sin6' T 3 I \ , n ! \ 

— o \R\}^aj R\J^b) “t“ Ron\l^a) 

Papp^ 


(3) 


Ron ()] • 


Since r oc At oc . Hence, the determination of 

time lags from the light curves at various frequencies can be used 
to determine the index kr- 

Q,ro is a measure of the core-position offset in pc(GHz)^/*’’' 
given by lLobanov|1998[( 


= 4.85 X 10 


— 9 AVnoasDL 

(1 + zf 


V 


1/kr 1/fcr 


1 / kr 


1/kr 


(4) 


where Armas = 's the core-position offset (in milli- 

arcseconds) between two frequencies v\ and V 2 (in GHz) in terms 
of the proper motion /r and Dr is the luminosity distance in pc. 
The distance of VLBI core from the base of the jet at frequency v 
is given by jLobanov|I998[l 


r core (r^) 


sin 6 


(5) 


Using Eqn. (43) from [Hirotani[ ( |2005l along with the assumption 
of equipartition between magnetic field energy density and particle 
energy density, the magnetic field strength in Gauss at 1 pc from 
the jet base is given by 


Hi 0.025 


f (1 + z f 

\ S'^i.psiv? 6 


1/4 


( 6 ) 


where 6 is the jet viewing angle, ip is the jet half-opening angle, 5 
is the Doppler factor, 2 is the cosmological redshift and Q.rv is the 
core-position offset defined in Eqn. Using the above equations 
and the assumption that m = 1 for an equipartition magnetic field. 


B 


core 


[v) = Bir 


-1 
core • 


(7) 
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Figure 1. Long-term variability light curves of 3C 454.3 in the 4.8 GHz - 
36.8 GHz frequency range. 


3 OBSERVATIONS AND ANALYSIS TECHNIQUE 

The long term light curves of 3C 454.3 at 4.8 GHz, 8.0 GHz and 
14.5 GHz are obtained from the University of Michigan Radio As¬ 
tronomical Observatory (UMRAO) ( [Aller et ar][1999[ (. The UM- 
RAO data at 4.8 GHz spans the period between 1978 and 2013, 
while the 8.0 GHz data spans from 1966 to 2013 and the 14.5 GHz 
data covers more than 35 years from 1973 to 2013. Radio mon¬ 
itoring at 22.2 and 36.8 GHz were carried out with the 22-m ra¬ 
dio telescope (RT-22) of the Crimean Astrophysical Observatory 
(CrAO) jVolvach|2006t . We use modulated radiometers in combi¬ 
nation with the registration regime “ON-ON” for collecting data 
from the telescope ( [Nesterov, Volvach & Strepka|2000[ l. Observa¬ 
tions at 37.0 GHz were made with the 14 m radio telescope of Aalto 
University Metsahovi Radio Observatory in Einland. Data obtained 
at Metsahovi and RT-22 were combined in a single array to supple¬ 
ment each other. A detailed description of the data reduction and 
analysis of Metsahovi data is given in [Teraesranta et ar]jl998[ l. 

The light curves of 3C 454.3 from 4.8 GHz to 36.8 GHz are 
displayed in Eig.[T] Long term monitoring is necessary for studying 
and analyzing various properties of AGNs such as the dynamical 
evolution, the inner sub-parsec structure, radiation mechanisms and 
location of radiating regions. 

The light curve is subjected to pre-processing in order to deter¬ 
mine the positions of the maxima and minima which represent the 
boundaries of the flaring portions. A Gaussian filtering with a suit¬ 
able parameter is applied to filter high frequency components and 
hence determine the extrema in the light curve. Then, a Gaussian 
function of the form 

( 8 ) 

is used to represent the shape of each flare region in the light curve. 
The fit is carried out from a base line at the minimum amplitude 
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of the light curve. Initial values Aq, Trio and cro are taken to be the 
maxima values, maxima positions and half the differences between 
consecutive minima (assuming that a flare would be symmetric) re¬ 
spectively. As Ao, Trio and cro are for a smoothed version of the 
raw light curve, trial values for A, m and a are generated based 
on regions around these initial values. These trial values account 
for possible values in a region around the flare peak amplitude, po¬ 
sitions and the differences in the original light curve. A matrix of 
{Ai,Tfii, (Ji) is constructed for all combinations of trial values of 
each peak. Each three parameter combination is fed into eqn. © 
along with the appropriate t spanning each flare. This procedure is 
cycled through the entire light curve such that all flares are covered 
by the Gaussian fit. For each Gaussian section covering a flare 

2 ^ [Xiitj) - y{Ai,Tni,cji)f 
^ VsiTmnce{xi{ti)) 

is determined for each section between the minima and the best 
fit for that section would be that combination {Ai,TTii,ai) which 
gives the least value of Xi ■ The 1 cr errors are determined for each 
parameter based on the table of Xi generated as the distance from 
the best fit minimum x^ value. 

The single Gaussian fit procedure is highly advantageous as it 
provides a good computational speed (few to tens of seconds de¬ 
pending on parameter grid size); the possibility of using a finely 
sampled grid of trial values for A, Tri and a with a typical grid size 
of 40000-50000 combinations; better fit values as evident from 
the data fits and the residuals; ease of procedure in determination 
of best fit parameters and associated errors. The light curves and 
the best fit Gaussian summed over each segment along with the fit 
residuals are presented in Figs. Hina and a The residual given 
by [y{Ai,TTii,(Ji) — a;i(fi)]/Standard Deviation(xi(ti)) is also 
calculated with the best fit {Ai,Tni, ar) for each Gaussian segment 
and plotted in the lower panels in the above figures. 

The light curves were split into two segments, one from the 
start of the observations to 2007.0 and the other from 2007.0 to 
the end of the observations. This was done to check for consis¬ 
tency of results between the first segment and the second segment, 
to determine any evolution in the fit parameters and due to the na¬ 
ture of the flaring in the light curve i.e. the flares before 2007.0 
are of lower amplitude compared to the post 2007.0 flares. If the 
entire light curve is run through the Gaussian fit procedure, it is 
likely to miss out prominent features in the first segment leading to 
incorrect results. We then conduct a time series analysis of each 
light curve segment using the Fourier periodogram analysis. As 
the periodogram is only applicable to evenly sampled light curves, 
we performed a linear interpolation and sampled the light curve 
at regular intervals of 0.1 days, similar to the procedure followed 
in, e.g. IGonzalez-Martln & Vaugh^j2012| l; |Mohan & Mangalam| 
( |2014^ . This evenly sampled light curve is used to determine the 
power spectral density (PSD) shape and any statistically significant 
quasi-periodic oscillations (QPOs). The normalized periodogram is 
given by (e.g. |Vaughan et al.|2003^ , 

o A+ 

Pifj) = (10) 

where At = 0.1 days is the sampling time step, x is the mean of 
the light curve x{tk) of length N points, and \F{fj)\ is its discrete 
Fourier transform evaluated at frequencies fj = j/{NAt) with 
j — 1, 2,.., (Af/2 —1). We fit the periodogram with two competing 
models, a power law and a bending power law, in order to describe 
the PSD shape. The power law model is given by 

i{fj) = Af; + c, ( 11 ) 


with amplitude A, slope s, and a constant Poisson noise C. This 
model fits optical/ultra-violet and X-ray data reasonably well as 
inferred from previous studies (e.g. [Gonzalez-Martin & Vaughan] 
|2012[ |Mohan & Mangalam|2014^ and could thus represent broad¬ 
band variability across a wide range of Fourier frequencies in mul¬ 
tiple wavelengths. The bending power law model is given by 

Hfj) = Mi" (1 + + c. ( 12 ) 


with amplitude A, slope s, bend frequency /(,, and a constant 
Poisson noise C. The model is a special case of the generalized 
bending knee model ( [McHardy et al.|2004^ . The maximum likeli¬ 
hood estimator (MLE) method is used to determine model param¬ 


eters in which the log-likelihood function (e.g. Emmanoulopoulos, 


McHmdy^P^adakisj20]JJjMohan, Mangalam & Chattopadhyay 
2014[ Mohan & Mangalam 2014| | defined by 


SiOk) = -2^(ln(I(/„0fe)) + P(/,)//(/„0O), (13) 

i=i 


is first determined. Here, I{fj,9k) are the power law or bend¬ 
ing power model with parameters 9k- Determining 9k which mini¬ 
mize S yields the maximum likelihood values. The log-likelihood 
S is determined for a large number of combinations of the pa¬ 
rameters 9k for a given model. The parameter combination which 
gives a unique global minimum Smin yields the best fit. For Si cor¬ 
responding to each unique combination of the parameters 9k, as 
AS = Si — S'min are approximately Xk distributed, the cumulative 
distribution function of the Xk distribution is used to estimate pa¬ 
rameter confidence intervals. The AS determined from the param¬ 
eter combinations within a desired confidence interval are grouped. 
Parameter ranges that they correspond to are used to determine con¬ 
fidence intervals of the 9k used in a given model. 

For model selection, we use the Akaike information criteria 
(AIC) (e.g. [Burnham & Anderson|2004| which measures the infor¬ 
mation lost when a model is fit to the data. The model with least 
information loss (least AIC) is the best fit. A likelihood function L 
which is proportional to the probability of the success of the model 
in describing the PSD shape is evaluated for each model. The AIC 
and likelihood are defined by 


AIC = S[9k) A ‘^Pk, (14) 

24/Gmin(model i) f^inin(iiull) 5 

L(model i|data) = 


where pk, the penalty term, is the number of parameters 9k used in 
the model, and L(model i|data) is the likelihood of model i given 
the data. Models with Ai ^ 2 can be considered close to the best 
fit, those with 4 A; ^ 7 are considerably less supported, and 
those with Ai > 10 cannot be supported (jBumham & Anderson| 


2004i. As the residual 7 (/j) = P{fj)/I{fj) follows the Xo distri¬ 


bution, the area under the tail of the probability density function of 
the X 2 distribution (gamma density r(l, 1 / 2 ) = exp (—a;/ 2 )/ 2 ) 
gives the probability e that the power deviates from the mean at 
a given frequency and is measured in units of standard deviation 
given by 7 ^. We must correct 7 ^ to account for the K number of 
trial frequencies at which the periodogram is evaluated and so is 
given by (|Vaughan|2005|l, 


7 e = -21n[l - (1 - e) ]. 


(15) 


Once e (e.g. 0.95, 0.99) is specified, 7 ^ is calculated and multiplied 
with l{fj) to give the significance level used to identify outliers in 
the periodogram that could indicate the presence of a QPO. 


































Frequency dependent core shifts and parameter estimation for the blazar 3C 454.3 5 



1985 1990 1995 2000 2005 

Time (Years) 



1985 1990 1995 2000 2005 

Time (Years) 



1980 1990 2000 

Time (Years) 


Figure 2. Segment 1 (beginning of observations to ~ 2007.0) of the 
4.8 GHz, 8.0 GHz and 14.5 GHz light curve flares fit with a piecewise 
Gaussian function. The residual in the lower panels are calculated as 
[y{Ai,Wii, cTi) — Xi(ii)]/Standard Deviation(xi(ti)). 
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Figure 3. Segment 1 (beginning of observations to ~ 2007.0) of the 22.0 
GHz and 36.8 GHz light curve flares fit with a piecewise Gaussian func¬ 
tion. The residual in the lower panels are calculated as [y(Ai,rni,ai) — 
Xi(ti )]/Standard Deviation(xi(t;)). 


4 RESULTS AND DISCUSSION 

For a flare X observed in multiple frequencies v, rfix,vi — 
rnx,v 2 ~ where Atx is the time lag between the flare in 

ui and ZZ 2 with vi > V 2 - The parameters derived from the Gaus¬ 
sian fit procedure discussed in section are presented in Table 
for the first segment and in Table|^for the second segment. In both 
tables, column 1 presents the flare nomenclature; column 2 gives 
the observed frequency, column 3 gives the maximum amplitude 
of the flare at a particular frequency; column 4 gives the epoch of 
the maximum flare amplitude; column 5 gives the full width at half 
maxima obtained from Gaussian fit for the particular flare; column 
6 gives the calculated time lags and column 7 gives the spectral 
index for each flare. 

We fit time 1 ms (Atx) derived from the Gaussian fits de¬ 
scribed in section M with a function At = + b using a 

weighted non-linear least squares method. The median of each pa¬ 
rameter and the associated median deviation are reported below. 
Those portions that were not fit well due to either too small a num¬ 
ber of points or unstable numerical results are excluded in the cal¬ 
culation of the parameters. The results of At versus v are plotted 
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Figure 4. Segment 2 (2007.0 to end of observations) of the 4.8 GHz, 8.0 
GHz and 14.5 GHz light cui've flares fit with a piecewise Gaussian func¬ 
tion. The residual in the lower panels are calculated as [y(Ai,mi, (7^) — 
iCi(fi)]/Standard Deviation(xi(ti)). 
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Figure 5. Segment 2 (2007.0 to end of observations) of the 22.0 GHz 
and 36.8 GHz light curve flares fit with a piecewise Gaussian func¬ 
tion. The residual in the lower panels are calculated as [y(Ai, (jj) — 
iCj(fi)]/Standard Deviation(xi(ti)). 


in Fig. for segment 1 and Fig. for segment 2 and the values of 
the fit parameters are: 

(i) First segment: beginning of observation - 2007.0 

b = -0.28 ± 0.19 t/r 
a = 5.51 ± 1.03 yr 
kr = 0.97 ±0.24; 

(ii) Second segment: 2007.0 - end of observation 

b = —0.06 ± 0.02 yr 
a = 3.03 ± 0.87 yr 
kr = 1.22 ±0.33. 

In both light curve segments, the inferred kr are similar and 
consistent with equipartition between particle energy density and 
magnetic field energy density. From the above fits, we have infor¬ 
mation on the amplitude and hence the peak flux for each flare at the 
observed frequencies. As S oc v°‘, the amplitude versus frequency 
data is fit with a function A = Al to determine the spectral index 
a for each flare. The fit is carried out only for flares with four data 
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Figure 6. Time lags At as a function of v for segment 1 light curves. A fit 
to At = + b yields b = —0.28 di 0.19 yr, a = 5.51 di 1.03 yr 

and kr = 0.97 di 0.24. The kr value indicates consistency 
of the equipartition between magnetic field energy density and the particle 
kinetic energy. 
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Figure 7. Time lags At as a function of p for segment 2 light curves. A fit 
to At = + b yields smaller magnitudes for b = —0.06 it 0.02 yr 

and a = 3.03 it 0.89 yr and a consistent = 1.22 it 0.33 

when compared with the fit to segment 1 light curves. 


points or more. We find a in the range —0.24 and 1.52 for segment 

1 and 1.16 and 1.31 for flares Q and S of segment 2. The calculated 
a are presented in Tahle[T]for segment 1 and in Tablej^for segment 
2. In a majority of cases (7/11), we infer a positive a indicating 
that there is in general an increase in amplitude with frequency for 
that flare. For two flares (C and K), we obtain the opposite trend, 
i.e. decreasing amplitude with increasing frequency, while for the 
remaining two (A and I) there is no significant dependence. 

As the equipartition condition is satisfied in both segments, 
we use eqns. and 0 to determine the magnetic field strength 
at 1 pc from the jet base and at the core respectively; B\, Score, 
the core position offset from eqn. 0 and the distance of the VLBI 
core from the base of the jet from eqn. 0 are calculated and pre¬ 
sented in Table for the first segment and Table for the sec¬ 
ond segment. Using 9 = 1.3°, cf> = 0.8°, <5 = 24.6, average 
proper motion ^ = 0.12 mas/yr ( [Jorstad et al. [2010) 1, 2 = 0.859 
and D = 5.489 Gpc (Ha = 0.73 for the vacuum energy den¬ 
sity, flm = 0.27 for the matter energy density and Hubble con¬ 
stant Ho = 71 km/s/Mpc), we obtain weighted mean values of 
Bi = 0.49 ± 0.18 G for segment 1, Bi = 0.36 ± 0.05 G for seg¬ 
ment 2, Score = 47 ± 16 mG for segment 1 and Score = 23 ± 7 
mG for segment 2. We obtain a weighted mean Q,ru = 6.59 ± 3.09 
pc GHz^/*^'' when averaged over all frequency pairs in the first seg¬ 
ment and Hri/ = 4.25 ± 0.78 pc in the second segment. 

The core radius rcore is expected to decrease with increasing i/. For 
9 out of 20 flares, namely, B, F, H, L, N, O, P, R, and T, no con¬ 
clusion can be drawn as the number of data points is two or one. 
The expected trend is inferred in 7 of the remaining 11 flares with 
the removal of one of the data points, the 4.8 GHz (flare K), 14.5 
GHz (flare E) and 22 GHz (flares A, C, G, I and S). The remain¬ 
ing 4 flares (D, J, M and Q) show the reverse trend. This indicates 
that the synchrotron opacity model may not be valid in all cases. 
In both tables, column 1 presents the flare nomenclature; column 

2 gives the observed frequency, column 3 gives the core position 
offset in mas; column 4 gives the core position offset in pcGHz; 
column 5 gives the distance of the emitting core from the base of 
the jet; column 6 gives the magnetic field strength in Gauss at 1 
pc and column 7 gives the magnetic field strength in Gauss at the 
emitting core. 

In the time series analysis of these light curve segments using 
the periodogram, the power law model describes the PSD shape 
better in 4/10 segments. The power law slope ranges between —1.6 
and —3.0 across all segments, consistent with that inferred for opti¬ 
cal and X-ray light curve studies (e.g.|Gonzalez-Martln & Vaughan] 
|2012[ [Mohan & Mangalam||2014^ . There is no clear trend of in¬ 
crease or decrease of the slope with observation frequency indi¬ 
cating that the origin of the flares are from multiple propagating 
shocks as also indicated from the study on flare amplitude depen¬ 
dence on observation frequency. The bending power law model 
describes the PSD shape better in the remaining 6/10 segments. 
The power law slope at lower temporal frequencies ranges between 
—2.5 and —3.5. The slopes for a majority of the second segment 
light curves are consistently steep at —3.5; since the bending power 
law is the best fit for all second segment light curves, it indicates 
that the second segment flares originate from regions of similar 
size. The timescale associated with the bend frequency ranges only 
between 0.51 yr and 0.66 yr for these second segment light curves, 
supporting the inference from the slopes with further information 
that the emitting region is at a similar distance from the base of the 
jet. The periodograms of each segment with the best fit PSD mod¬ 
els are presented in Figs.[ 8 ][^| 10 |and[TT]and detailed results from 
the time series analysis are presented in Table[^ 
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3C 454.3 4.8 GHz: Segment 1 



3C 454.3 8.0 GHz: Segment 1 



3C 454.3 14.5 GHz: Segment 1 



Figure 8. Periodogram analysis of segment 1 (beginning of observations to 
~ 2007.0) of the 4.8 GHz, 8.0 GHz and 14.5 GHz light curve flares. The 
best fit model is the solid curve, the dashed curve above it is the 99 % sig¬ 
nificance contour which can identify statistically significant quasi-periodic 
components, the dot-dashed horizontal line is the white noise level and the 
plot below each periodogram panel shows the fit residuals. 


Flare 

Frequency 

Amplitude 

Position 

Width 

Time lag 

Spectral index 


(GHz) 

A(Jy) 

m (year) 

a (years) 

At (years) 

a 

A 

36.8 

9.28± 0.40 

1981.58±0.08 

0.86±0.19 

0.00 

-0.09±0.13 


22.0 

10.97±0.18 

1982.02±0.14 

0.82±0.39 

0.44±0.16 



14.5 

12.63±0.06 

1981.80±0.10 

0.96±0.51 

0.22±0.17 



8.0 

12.38±0.30 

1982.25±0.09 

0.96±0.36 

0.67 ±0.13 



4.8 

9.01±0.54 

1982.33±0.17 

1.02±0.15 

0.75±0.19 


B 

8.0 

2.90±0.46 

1984.74±0.14 

0.39±0.06 

0.00 



4.8 

3.39±0.10 

1985.09±0.16 

0.55±0.11 

0.35±0.21 


C 

36.8 

4.66±0.36 

1987.21±0.02 

0.95±0.33 

0.00 

-0.19 ±0.13 


22.0 

5.46±0.28 

1987.09±0.13 

1.54±0.43 

-0.12±0.13 



14.5 

7.16±0.02 

1987.33±0.06 

1.03±0.50 

0.12±0.14 



8.0 

7.62±0.42 

1987.43±0.03 

1.17±0.25 

0.22±0.07 



4.8 

6.98±0.46 

1987.76±0.06 

1.33±0.27 

0.55±0.07 


D 

36.8 

7.22±0.28 

1989.05±0.02 

0.45±0.27 

0.00 

0.12 ±0.02 


22.0 

6.00±0.38 

1989.30±0.18 

0.48±0.{)4 

0.25 ±0.18 



14.5 

6.37±0.42 

1989.36±0.13 

0.55±0.16 

0.31 ±0.22 



8.0 

6.21±0.52 

1989.40±0.03 

0.59±0.18 

0.35±0.13 



4.8 

5.47±0.14 

1989.10±0.01 

1.09±0.30 

0.05 ±0.03 


E 

36.8 

6.26±0.56 

1991.01±0.09 

0.23±0.12 

0.00 

1.52 ±0.23 


22.0 

2.45±0.24 

1991.16±0.14 

0.31±0.08 

0.16±0.17 



14.5 

1.67±0.32 

1991.37±0.15 

0.48±0.10 

0.37±0.21 



8.0 

0.86±0.28 

1991.44±0.03 

0.13±0.03 

0.44 ±0.15 


F 

36.0 

8.81±0.08 

1993.26±0.05 

0.36±0.15 

0.00 



22.0 

9.2±0.02 

1993.44±0.13 

0.98±0.31 

0.18±0.14 


G 

36.8 

12.41±0.01 

1994.05±0.15 

0.41±0.14 

0.00 

0.16 ±0.01 


22.0 

12.28±0.30 

1994.29±0.08 

0.63±0.17 

0.24±0.17 



14.5 

15.61±0.24 

1994.18±0.02 

0.92±0.25 

0.13±0.08 



8.0 

13.37±0.30 

1994.42±0.16 

0.77±0.27 

0.37±0.16 



4.8 

8.94±0.01 

1994.82±0.11 

1.06±0.37 

0.77±0.19 


H 

36.8 

7.50±0.10 

1995.26±0.18 

0.50±0.01 

0.00 



22.0 

8.01±0.42 

1995.42±0.07 

0.42±0.07 

0.16±0.19 



14.5 

7.87±0.20 

1995.57±0.19 

0.37±0.05 

0.31 ±0.20 


I 

36.8 

4.43±0.34 

1998.69±0.10 

0.48±0.17 

0.00 

-0.04 ± 0.09 


22.0 

6.10±0.46 

1998.57±0.09 

0.36±0.09 

-0.12±0.13 



14.5 

6.14±0.26 

1998.65±0.16 

0.41±0.06 

-0.04±0.18 



8.0 

5.60±0.14 

1998.73±0.17 

1.29±0.34 

0.04±0.23 



4.8 

4.79±0.52 

1998.95±0.06 

1.22±0.38 

0.26±0.18 


J 

36.8 

4.95±0.12 

2000.04±0.14 

0.77±0.15 

0.00 

0.08 ± 0.06 


22.0 

5.10±0.08 

2000.36±0.05 

0.31±0.13 

0.32±0.15 



14.5 

4.60±0.12 

2000.37±0.13 

0.91±0.27 

0.33±0.14 



8.0 

3.61±0.58 

2000.61±0.06 

0.28±0.10 

0.57±0.14 



4.8 

4.41±0.38 

2001.22±0.09 

0.78±0.26 

0.18±0.11 


K 

36.8 

4.37±0.50 

2001.36±0.12 

1.01±0.18 

0.00 

-0.24 ±0.17 


22.0 

4.61±0.04 

2001.38±0.14 

0.65±0.11 

0.02±0.18 



14.5 

5.52±0.10 

2001.56±0.03 

0.49±0.16 

0.20±0.14 



8.0 

4.20±0.36 

2002.29±0.20 

1.46±0.37 

0.93±0.20 



4.8 

4.21±0.58 

2002.85±0.01 

1.23±0.31 

1.49±0.20 


L 

36.8 

5.27±0.10 

2003.20±0.09 

0.35±0.08 

0.00 



22.0 

4.38±0.20 

2003.27±0.12 

0.64±0.16 

0.07 ±0.15 



14.5 

5.05±0.36 

2003.24±0.01 

0.77±0.17 

0.04 ±0.12 


M 

36.8 

17.55±0.02 

2006.06±0.03 

0.26±0.37 

0.00 

1.05 ±0.06 


22.0 

8.07±0.44 

2006.24±0.18 

0.34±0.03 

0.16±0.18 



14.5 

5.08±0.50 

2006.35±0.05 

1.09±0.15 

0.25 ±0.19 



8.0 

4.38±0.22 

2006.20±0.01 

0.58±0.45 

0.14±0.05 



Table 1. Gaussian fit based parameters, time lags and spectral index for 
segment 1 light curves. 


5 CONCLUSIONS 

The core shift effect in the parsec scale jet of the blazar 3C 454.3 
was studied in the 4.8 GHz - 36.8 GHz radio wavelengths. The 
inferred time delays Af for flares in these frequencies are in agree¬ 
ment with the frequency dependence of synchrotron emission from 
the core (surface with optical depth t — 1) jKonigl|1981| l. From 
the fit At oc where u are the observation frequencies, we 

obtain a weighted mean kr = 1.10 ± 0.18, consistent with the 
equipartition between the magnetic field energy density and the 
particle energy density. The large scale magnetic field thus domi¬ 
nates the kinematics of parsec scale jets as its decay is slower com¬ 
pared to the decay in the particle number density as Score oc 
and N oc r~^. From the fit to the flare amplitudes A oc 
we infer spectral indices a in the range —0.24 and 1.52 for both 
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Flare 

Frequency 

Amplitude 

Position 

Width 

Time lag 

Spectral index 


(GHz) 

^(Jy) 

fn (year) 

cr (years) 

At (years) 

a 

N 

14.5 

1.08±0.12 

2006.85±0.08 

0.28±0.12 

0.00 



8.0 

2.25±0.54 

2007.19±0.01 

0.23±0.06 

0.34±0.08 



4.8 

2.15±0.56 

2007.29±0.13 

0.50±0.15 

0.44±0.13 


O 

36.0 

11.19±0.22 

2007.66±0.01 

0.15±0.05 

0.00 



22.0 

9.50±0,28 

2{)08.06±0.12 

0.23±0.02 

0.40±0.12 



14.5 

3.25±0.36 

2{)08.18±0.14 

0.41±0.20 

0.52±0.14 


P 

36.0 

11.61±0.46 

2{)07.96±0.06 

0.17±0.03 

0.00 



4.8 

1.98±0.24 

2008.61±0.23 

0.48±0.04 

0.65±0.24 


Q 

36.0 

24.44± 0.60 

2008.68±0.07 

0.26±0.09 

0.00 

1.16 ±0.06 


22.0 

15.36±0.06 

2008.92±0.04 

0.27±0.26 

0.24±0.08 



14.5 

9.31±0.18 

2008.84±0.04 

0.33±0.14 

0.16±0.04 



8.0 

4.38±0,08 

2008.90±0.16 

0.60±0.21 

0.22±0.16 



4.8 

3.57±0.54 

2009.31±0.09 

1.02±0.15 

0.53±0.10 


R 

36.0 

32.43 ±0.24 

2010.31±0.02 

0.34±0.09 

0.00 



22.0 

31.85±0,34 

2010.54±0.01 

0.36±0.12 

0.23±0.02 



14.5 

22.34±0.06 

2010.63±0.07 

0.34±0.15 

0.32±0.07 


S 

36.0 

47.99±0.12 

2010.96±0.01 

0.28±0.27 

0.00 

1.31 ±0.24 


22.0 

35.79±0.16 

2010.91±0.01 

0.13±0.06 

-0.05 ±0.01 



14.5 

12.63±0.06 

2011.03±0.02 

0.96±0.51 

0.07±0.02 



8.0 

14.27±0.24 

2011.22±0.32 

0.75±0.03 

0.26±0.32 



4.8 

7.80±0.28 

2011.46±0.15 

0.85±0.27 

0.50±0.35 


T 

36.0 

9.64± 0.04 

2013.78±0.01 

0.50±0.49 

0.00 



22.0 

3.39±0.01 

2013.81±0.12 

0.30±0.92 

-0.03±0.12 



Table 2. Gaussian fit based parameters, time lags and spectral index for 
segment 2 light curves. 


segments. As the shock propagates downstream, it loses energy 
through interactions with the plasma. Thus, the amplitude of flares 
at lower frequencies is expected to be reduced. The opposite trend, 
i.e., decreasing amplitude with increasing frequency, can occur if 
the core is excited by shocks that have sufficient energy to cause 
flaring at lower frequencies but are not energetic enough to cause 
excitation upstream at higher frequencies. In our study, we obtain 
both positive and negative a. The flaring activity can thus be at¬ 
tributed to multiple propagating shocks. 

We infer a weighted mean Bi = 0.5 ± 0.2 G and Score = 
46± 16 mG considering both segments. These results are consistent 
within error bars of Bi = 0.493 G obtained in |Chai, Cao & Gu] 
( [2012^ and Score = 0.04±0.02 G obtained in |Kutkin et al. ( |2014| >. 
We obtain a weighted mean Qri^ ~ 6.4 ± 2.8 pc when 

averaged over all frequency pairs in both segments. The source was 
also studied by |Pushkarev et al^j2012^ where estimates of fin/, Si 
and Score for the 15/8 GHz frequency pair are 22 pcGHz (kr = 1 
was assumed), 1.13 G and 0.06 G respectively. The inferred Si and 
Score range in our study are consistent with these earlier estimates. 
The shift /r = 0.70 mas/yr was assumed in |Kutkin et al.|j2014^ 
and as a result, t}rv = 43 ± 10 pc GHz was obtained. As it was 
noted that the chosen value is 2 — 8 times higher than that obtained 
in previous studies, we used a lower limit of ^ — 0.12 mas/yr of 
a range 0.12 - 0.53 m as/yr jj orstad e t al.||2005] l at 43 GHz, and 
0.3 mas/yr at 15 GHz ( [Lister et al, 12013) to estimate which 
is the reason for our estimate being ~ 4 — 15 times smaller than 
that obtained previously. Based on the statistical trend shown by 
the rcore as a function of v, we infer that the synchrotron opacity 
model may not be valid for all cases. 

From the time series analysis we obtain typical power law 
slopes in the range —1.6 to —3.5 for the power law as well as bend¬ 
ing power law PSD shapes at lower temporal frequencies. The anal¬ 
ysis indicates that the first segment light curves are consistent with 
multiple shock excitation events. The bending power law model is 
inferred to be a better PSD fit with a —3.5 slope and consistent bend 
timescales ranging between 0.51 yr and 0.66 yr for the second seg- 


Flare 

Frequency 

Ar 


Tcore 

Bfpc 

Bcore 


(GHz) 

(^as) 

(pc GHzl/*:^) 

(pc) 

(G) 

(mG) 

A 

4.8 

90±16 

3.97±0.74 

34.84±0.S3 

0.37±0.05 

11±2 


8.0 

80±20 

6.63±1.67 

34.47±1.82 

0.55±0.10 

16±4 


14.5 

26±19 

5.16±3.72 

14.54±4.02 

0.45±0.25 

31±23 


22.0 

53±10 

23.76±4.52 

43.62±4.90 

1.43±0.20 

33±6 

B 

4.8 

42±17 

1.85±0.76 

16.26±0.83 

0.21±0.06 

13±5 

C 

4.8 

66±8 

2.91±0.39 

25.55±0.45 

0.30±0.03 

12±2 


8.0 

26±17 

2.18±1.41 

11.32±1.52 

0.24±0.12 

21±14 


14.5 

14±16 

2.81±3.13 

7.93±3.38 

0.29±0.24 

36±41 


22.0 

14±2 

6.48±0.91 

11.90±0.98 

0.54±0.06 

45±6 

D 

4.8 

6±16 

0.26±0.71 

2.32±0.76 

0.05±0.10 

21±56 


8.0 

42±26 

3.47±2.15 

18.00±2.33 

0.34±0.16 

19±12 


14.5 

37±22 

7.27±4.30 

20.49±4.66 

0.59±0.26 

29±17 


22.0 

30±2 

13.50±0.93 

24.79±1.02 

0.94±0.05 

38±3 

E 

8.0 

53±25 

4.36±2.07 

22.63±2.24 

0.40±0.14 

18±9 


14.5 

44±20 

8.68±3.91 

24.46±4.24 

0.67±0.23 

27±13 


22.0 

19±11 

8.64±4.95 

15.86±5.36 

0.67±0.28 

42±25 

F 

22.0 

22±6 

9.72±2.71 

17.85±2.93 

0.73±0.02 

41± 12 

G 

4.8 

92±19 

4.07±0.87 

35.77±0.97 

0.38±0.06 

11±2 


8.0 

44±10 

3.67±0.84 

19.03±0.91 

0.35±0.06 

19±4 


14.5 

16±20 

3.05±3.91 

8.59±4.23 

0.31±0.29 

36±47 


22.0 

29±18 

12.96±8.10 

23.80±8.77 

0.91 ±0.43 

38±25 

H 

14.5 

37±23 

7.27±4.50 

20.49±4.87 

0.59±0.27 

29±18 


22.0 

19±22 

8.64±9.90 

15.86±10.71 

0.67±0.57 

42±53 

I 

4.8 

31±28 

1.38±1.24 

12.08±1.34 

0.17±0.11 

14±13 


8.0 

5±22 

0.40±1.82 

2.06±1.96 

0.07±0.23 

32±150 


14.5 

5±16 

0.94±3.13 

2.64±3.38 

0.13±0.32 

48±160 


22.0 

14±12 

6.48±5.40 

11.90±5.84 

0.54±0.34 

45±38 

J 

4.8 

22±17 

0.95±0.75 

8.36±0.81 

0.13±0.08 

15±12 


8.0 

68±17 

5.64±1.42 

29.32±1.55 

0.49±0.09 

17±4 


14.5 

40±18 

7.74±3.52 

21.81±3.81 

0.62±0.21 

28±13 


22.0 

38±17 

17.28±7.66 

31.73±8.28 

1.13±0.37 

36±16 

K 

4.8 

179±24 

7.88±1.15 

69.21±1.32 

0.62±0.07 

9±1 


8.0 

112±17 

9.21±1.44 

47.84±1.59 

0.70±0.08 

15±2 


14.5 

24±22 

4.69±4.30 

13.22±4.65 

0.42±0.29 

32±30 


22.0 

2±14 

1.08±6.30 

1.98±6.81 

0.14±0.62 

71±456 

L 

14.5 

5±18 

0.94±3.51 

2.64±3.81 

0.13±0.36 

48±186 


22.0 

8±11 

3.78±4.95 

6.94±5.36 

0.36±0.35 

52±70 

M 

8.0 

17±23 

1.39±1.90 

7.20±2.05 

0.17±0.17 

24±33 


14.5 

30±22 

5.86±4.30 

16.52±4.65 

0.50±0.28 

30±22 


22.0 

19±4 

8.64±1.81 

15.86±1.96 

0.67±0.10 

42±9 


Table 3. Core position offsets, distance from jet base and magnetic field 
strengths inferred for segment 1 light curves. 


Flare 

Frequency 

Ar 

flru 

Tcore 

B\pc 

Bcore 


(GHz) 

(mas) 

(pc GHz’-'^^’') 

(pc) 

(G) 

(G) 

N 

4.8 

53±10 

1.80±0.36 

22.11±0.87 

0.21±0.03 

9±2 


8.0 

41±10 

2.40±0.60 

19.44±1.40 

0.26±0.05 

13±3 

0 

14.5 

62±6 

8.01±0.80 

39.77±1,90 

0.63±0.05 

16±2 


22.0 

11±1 

3.03±0.29 

10.69±0.67 

0.31±0.02 

29±3 

P 

4.8 

78±7 

2.67±0.30 

32.66±0.7S 

0.28±0.02 

8±1 

Q 

4.8 

76±19 

2.59±0.67 

31.66±1.60 

0.27±0.05 

9±2 


8.0 

26±5 

1.56±0.30 

12.58±0.71 

0.19±0.03 

15±3 


14.5 

19±10 

2.47±1.29 

12.24±2.9S 

0.26±0.10 

21±11 


22.0 

28±8 

7.73±2.25 

27.31±5.22 

0.62±0.13 

23±7 

R 

14.5 

3S±2 

4.93±0.29 

24.47±0.71 

0.44±0.02 

18±1 


22.0 

29±2 

8.07±0.59 

28.50±1.38 

0.64±0.03 

22±2 

S 

4.8 

60±38 

2.05±1.31 

25.13±3.04 

0.23±0.11 

9±6 


8.0 

31±2 

1.84±0.14 

14.86±0.36 

0.21±0.01 

14±1 


14.5 

8±1 

1.08±0.13 

5.35±0.31 

0.14±0.01 

26±3 


22.0 

6±1 

1.68±0.28 

5.94±0.66 

0.20±0.02 

33±6 

T 

22.0 

4±1 

1.01±0.28 

3.56±0.65 

0.13±0.03 

38±10 


Table 4. Core position offsets, distance from jet base and magnetic field 
strengths inferred for segment 1 light curves. 
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3C 454.3 22.0 GHz: Segment 1 


3C454.3 4.8 GHz: Segment 2 



3C 454.3 36.8 GHz: Segment 1 



Figure 9. Periodogram analysis of segment 1 (beginning of observations to 
~ 2007.0) of the 22.0 GHz and 36.8 GHz light curve flares. The best fit 
model is the solid curve, the dashed curve above it is the 99 % significance 
contour which can identify statistically significant quasi-periodic compo¬ 
nents, the dot-dashed horizontal line is the white noise level and the plot 
below each periodogram panel shows the fit residuals. 


ment light curves indicating that the emitting cores for these flares 
are similar in size and distance from the jet base. 

Important objectives met in this study include the demonstra¬ 
tion of the computational efficiency and statistical basis of the anal¬ 
ysis technique involving the piecewise Gaussian fit, which is auto¬ 
mated in the determination of relevant flares; the ability to use this 
analysis across a wide range of light curves; the check for and con¬ 
firmation of consistency of the results obtained for both light curve 
segments; the check for consistency with previously reported re¬ 
sults based on a single flare in |Kutkin et al.|p014^ and the evidence 
it provides for the core shift dependence on the frequency and its 
usefullness in the determination of a variety of jet diagnostics in the 
region close to the resolving limit of VLBI (in the near vicinity of 
the central supermassive black hole). 




3C 454.3 14.5 GHz; Segment 2 



Figure 10. Periodogram analysis of segment 2 (2007.0 to end of observa¬ 
tions) of the 4.8 GHz, 8.0 GHz and 14.5 GHz light curve flares. The best fit 
model is the solid curve, the dashed curve above it is the 99 % significance 
contour which can identify statistically significant quasi-periodic compo¬ 
nents, the dot-dashed horizontal line is the white noise level and the plot 
below each periodogram panel shows the fit residuals. 
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3C 454.3 22.0 GHz: Segment 2 



3C 454.3 36.8 GHz: Segment 2 



Figure 11. Periodogram analysis of segment 2 (2007.0 to end of observa¬ 
tions) of the 22.0 GHz and 36.8 GHz light curve flares. The best ht model is 
the solid curve, the dashed curve above it is the 99 % signihcance contour 
which can identify statistically signihcant quasi-periodic components, the 
dot-dashed horizontal line is the white noise level and the plot below each 
periodogram panel shows the ht residuals. 
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Observation 

PSD 

PSD Fit parameters 


AIC 

Model 

Frequency & 

mode! 





likelihood 

Segment 


log(A) 

Q 

log(fb) 



4.8 GHz; Seg. 1 

PL 

-3.3 ±0.1 

-2.7 ± 0.2 


1946.55 

1.00 


BPL 

-1.7 ±0.2 

-2.9 ± 0.3 

-0,88 ±0.21 

1947.78 

0.54 

4.8 GHz; Seg. 2 

PL 

-2.8 ± 0.2 

-2.6 ± 0.4 


339.74 

0.05 


BPL 

-2.0 ± 0.2 

-3.5 ± 0.4 

-0.19 ±0.09 

333.72 

1.00 

8.0 GHz; Seg. 1 

PL 

-2.9 ±0.1 

-2.5 ± 0.2 


4301.24 

0.17 


BPL 

-1.6 ±0.1 

-2.5 ± 0.3 

-0.88 ±0.16 

4297.69 

1.00 

8.0 GHz; Seg. 2 

PL 

-2.6 ± 0.2 

-3.0 ± 0.3 


795.54 

10-^ 


BPL 

-1.7 ±0.2 

-3.5 ± 0.2 

-0.20 ±0.01 

782.67 

1.00 

14.5 GHz; Seg. 1 

PL 

-2.5 ±0.1 

-2.3 ± 0.2 


1937.31 

1.00 


BPL 

-1.6 ±0.1 

-3.5 ± 0.3 

-0.32 ±0.12 

1944.49 

0.03 

14.5 GHz; Seg. 2 

PL 

-2.1 ±0.2 

-3.0 ± 0.3 


294.36 

0.01 


BPL 

-1.3 ±0.1 

-3.5 ± 0.2 

-0.20 ±0.01 

285.16 

1.00 

22.0 GHz; Seg. 1 

PL 

-2.4 ±0.1 

-2.4 ± 0.3 


1528.90 

1.00 


BPL 

-1.7 ±0.2 

-3.5 ± 0.2 

-0,23 ±0.12 

1544.59 

10-“ 

22.0 GHz; Seg. 2 

PL 

-1.6 ±0.1 

-2.5 ± 0.3 


754.88 

0.02 


BPL 

-0.8 ± 0.2 

-3.0 ± 0.3 

-0.18 ±0.01 

746.90 

1.00 

36.8 GHz; Seg. 1 

PL 

-2.0 ±0.1 

-1.9 ±0.2 


1242.52 

1.00 


BPL 

-1.5 ±0.2 

-3.3 ± 0.4 

-0.09 ±0.17 

1249.49 

0.03 

36.8 GHz; Seg. 2 

PL 

-1.5 ±0.1 

-2.9 ± 0.3 


270.35 

0.38 


BPL 

-0.8 ±0.1 

-3.3 ± 0.3 

-0,29 ±0.01 

268.41 

1.00 


Table 5. Results from the parametric PSD models fit to the periodogram. 
Columns 1-7 give the observation frequency and segment, the model (PL: 
power law -H constant noise, BPL: bending power law -H constant noise), the 
best-fit parameters log(V), slope a and the bend frequency fi, with their 
95% errors derived from AS, the AIC and the likelihood of a paificular 
model. The best fit PSD is highlighted. 
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